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ABSTRACT 



Context. Attenuation of the TeV y-ray flux from distant blazars through pair production with extragalactic background light leads 
to the development of electromagnetic cascades and subsequent, lower energy, GeV secondary y-ray emission. Due to the deflection 
of VHE cascade electrons by extragalactic magnetic fields (EGMF), the spectral shape of this arriving cascade y-ray emission is 
dependent on the strength of the EGMF. Thus, the spectral shape of the GeV-TeV emission from blazars has the potential to probe the 
EGMF strength along the line of sight to the object. Constraints on EGMF previously derived from the gamma-ray data suffer from 
an uncertainty related to the non-simultaneity of GeV and TeV band observations. 

Aims. We investigate constraints on the EGMF derived from observations of blazars for which TeV observations simultaneous with 
those by Fermi telescope were reported. We study the dependence of the EGMF bound on the hidden assumptions it rests upon. 
Metliods. We select blazar objects for which simultaneous Fermi/LAT GeV and Veritas, MAGIC or HESS TeV emission have been 
published. We model the development of electromagnetic cascades along the gamma-ray beams from these sources using Monte 
Carlo simulations, including the calculation of the temporal delay incurred by cascade photons, relative to the light propagation time 
of direct y-rays from the source. 

Results. Constraints on EGMF could be derived from the simultaneous GeV-TeV data on the blazars RGB J0710+591, 
lES 0229-1-200, and lES 1218-1-304. The measured source flux level in the GeV band is lower than the flux of the expected cascade 
component calculated under the assumption of zero EGMF. Assuming that the reason for the suppression of the cascade component is 
the extended nature of the cascade emission, we find that B > 10"'^ G (assuming EGMF correlation length of > 1 Mpc) is consistent 
with the data. Alternatively, the assumption that the suppression of the cascade emission is caused by the time delay of the cascade 
photons the data are consistent with B > 10"" G for the same correlation length. 

Key words. Gamma rays: galaxies - Galaxies: active - BL Lacertae objects: general 



1. Introduction 

The presence of magnetic fields in galaxies and galaxy clus- 
ters plays a key role in present day astrophysical studies. 
Hovyever, t he ori g in of these fields rema i ns largely uncertain (se e 
Kronberg|([l994l) : lGrasso & RubinsteinI (l200ll) : IWidro"wl(l2002h : 



BecklT2()09l) for reviews). A commonly accepted hypothesis is 



that relatively strong galactic and cluster magnetic fields re- 
sult from the amplification of much weaker pre-existing "seed" 
fields via compression and turbulence/dyna mo amplification in 
the co urse of structure formation processes dKulsrud & Zweibeil 
120081) . 

The origin of these seed magnetic fields is unknown. 
It is possible that the seed fields are produced locally in 
(proto)galaxi es via th e so-called " Biermann ba ttery" mech- 
anism TPudri tz & SiiS 1989; Gnedi n et al.l l2000h . Otherwise, 
the seed fields might be of primordial origin, i.e. produced 
at the moments of phase transitions in the Early Universe 
dOrasso & Rubinsteinll200lt IWidrowll2002l) . Constraints on the 
nature of the seed fields could potentially be derived from 
the measurements of weak magnetic fields in the intergalactic 
medium which are not amplified by the action of different types 
of dynamos. 

The measurement of extremely weak magnetic fields in 
the voids of the Large Scale Structure (LSS) is a challenging 
task and up to now only upper bounds have been derived us- 



ing various techniques. The tightest upper bounds come from 
the search for the Faraday r otation of p olarizati on of radio 
emiss i on from distant quasars jKronberg & Simar d-NormandinI 
119761: iKronberg & Perrvlll982HBlasi et al.iri999l) and from the 
effect of magnetic field s on the anisotropy of Cosmic Microwa ve 
Background radiation jBarrow et alJll997l:lDu"rrer et alj|2000h . 

A new handle on the EGMF measure, using the cascade 
emission from blazars, is now emerging as an alterna- 
tive probe. In this method, multi-TeV y-rays from distant 
(> 100 Mpc) blazars attenuate through pair production 
interactions on the extragalactic background light (EBL), 
leading to the develo p ment o f elec tr omagnetic cascades 



(lAharonian e t alj [19941: iPlagal 119951: ICoppi & AharonianI 



1996: Nero nov & Semikozl [20071 Id'Avezac et al.l [2007 
Murase et al. 2008: Eunawanichavapant & Aharonian 
.Neronov & Semikoai2009;.Elviv et al...2009; .Dolag et al. 



20091: 
i,2009i) . 

The angular pattern of the secondary cascade emission from 
e^e" pairs deposited in the intergalactic medium through 
pair production interactions depends on the EGMF strength. 
The detection (non-detection) of the cascade emission signal 
from known TeV y-ray emitting blazars could result in the 
measurement of (lower bound on) the strength of the magnetic 
field in intergalactic space along the line of sight toward these 
blazars. The first application of this method f or deriving lower 
bounds on the EGMF have been carried out ( Neronov & Vovkl 
120101: lTavecchidl201flt [Polag et al.l[20T(it iDermer et al.ll2010h . 
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suggesting that a measure of the EGMF may finally soon be 
within reach. 

In the simplest settings, the lower bounds on EGMF at the 
level of 10 '^ to 10~'^ G (iNeronov & Vovkll201(]l: iTavecchiol 
I2OIOI: [Polag et al.]|2010h (depending on assumptions about in- 
trinsic blazar spectra) were derived. These bounds adopt a sim- 
plifying assumption that the measured blazar fluxes provide cor- 
rect estimates of the time-averaged fluxes in the GeV and TeV 
energy bands. This assumption could, in principle, be partially 
verified via a systematic monitoring of the sources simultane- 
ously in the TeV (for primary source emission) and GeV (for the 
cascade emission) energy bands. 

In this work, we consider a set of blazars observed simultane- 
ously both by Fermi/LAT and HESS or Veritas, in order to search 
for the cascade component of the GeV- TeV spectra of these 
sources. We find that in several cases, namely for the blazars 
RGB J0710+591 and lES 12184-304, the measured source flux 
in the GeV band is lower than the expected (minimal possi- 
ble) flux produced by the gamma-ray cascade in intergalactic 
space, calculated assuming zero magnetic field along the line of 
sight. This imposes a lower bound on the strength of magnetic 
field in the intergalactic medi um, similar to the bounds found for 
the cases of lES 02294- 200 dNeronov & Vovk||2010l:lTavecchio 
20To|: iDolag et al.ll20Tol) and lES 0347-121 dNeronov & Vovk 



2010D. However, contrary to the bounds from lES 02294-200 
and lES 0347-121, no additional assumptions about long-term 
stability of the source in the GeV and TeV bands are needed, 
because of the truly simultaneous nature of multi-band observa- 
tionfl 

Availability of the simultaneous data allows a study of two 
alternative possible reasons for the suppression of the cascade 
signal: dilution of the cascade flux due to the time delay of the 
cascade signal, following a period of enhanced activity of the 
source vs. suppression of the cascade signal contribution to the 
point source flux due to the extended nature of the cascade emis- 
sion. We show that adopting an (possibly extreme) assumption 
that the source gamma-ray activity periods of the sources are as 
short as ~ 2 yr (the time scale of Fermi operation), the lower 
bound on the magnetic field relaxes to ~ 10""^ G for the cases 
of RGB J0710+591 and lES 12184-304. 

When this paper was almost ready for publication a re- 
port on investigation of constraints on EGMF derived assum- 
ing suppression of the casca de emission due to the time delay 
of the cascade flux app eared ( iDermer et al. I I2OIOI) . The analysis 
of iDermer et alJ (1201 Ol) reUes on an unpubU shed VERITAS ob- 
servation of lES 02294-200 (|Perkinsll2010h which reveals that 
the flux and spectrum of the source did not chan ge on a 4-year 
l .ime span betw een HESS (lAharonian et al.ll2007h and VERITAS 
(|Perkinsll2010l) observations. Relying on this statement, we in- 
clude lES 02294-200 in our an alysis, to verify the results ob- 
tained by IDermer et al.l (l2010l) . We find that simplifying as- 
sumptions adopt ed in the ana l ytical modeling of electromag- 
netic cascade by IDermer et al.l (1201 Ol) led to a large underesti- 
mate of the lower b ound on EGMF. Correction of the result of 
IDermer et alj ( 1201 Ol) found using full Monte-Carlo simulation of 
the electromagnetic cascade we find a lower bound B > 10 G 
from the minimal possible time delay of the cascade signal in 
lES 02294-200. 



2. Selection of sources and data analysis 

Detailed calculation of the spectrum of GeV-band emission from 
electromagnetic cascade initiated by absorption of very-high- 
energy y-rays on the EBL requires the knowledge of the initial 
(unabsorbed) source flux in the TeV energy band. If the mag- 
netic field in the intergalactic medium is close to zero, the cas- 
cade GeV y-rays arrive almost simultaneously with the primary 
absorbed TeV y-rays, with only a small (~ 10 hr scale) mag- 
netic field-independent time delay related to the angular scat- 
ter of electrons and positr ons in the pair production process 
eronov & Semiko 3 I2OO9I) . This means that prediction of the 
level of the GeV-band cascade emission at a given moment of 
time requires the knowledge of the simultaneous TeV flux of the 
source. 

The Large Area Telescope (LAT) on board of Fermi satel- 
lite has performed continuous monitoring of the entire sky in 
the GeV band since August 4, 2008. This implies that any 
blazar observations in the TeV band, performed with the ground- 
based gamma-ray telescopes HESS, MAGIC or Veritas after 
August 4, 2008 automatically have simultaneous observations 
in the GeV band. Several blazar observations performed af - 
ter August 4, 2008 were reported: Mrk 501 ( lHuang|j2009|) . 
PKS 2155-304 ( lAharonian et al. 2009), 3C 66A at ( lRevesl2009l) . 
RGB J07 104-591 (lAcciari et al. 201 Ota) and_JES_m8±304 
dAcciariet al.| l2010cl) lES 0229-H 200 dDermer et al.1 120101: 
lPerkinsll201cil) and PKS 14244-240 (lAcciari et alj|2010al) . Our 
analysis of constraints on EGMF is based on the analysis of 
GeV- TeV band spectral properties of the sources listed above. 

For each of the sources with simultaneous GeV- TeV data, 
we derive the spectral characteristics in the 0.1-100 GeV band 
from FermijhPil data and combine them with the reported TeV 
band spectrum to produce the broad band spectra. In our anal- 
ysis of Fermi/LAI data we use the publicly available data from 
the Fermi/LAT from August 4, 2008 tifl November 30, 2010. 
We process the data using Fermi Science Tools package of the 
version v9rl7p(3. For the spectral extraction we use the un- 
binned likelihood analysis method, taki ng into account a ll the 
sources from the first year Fermi catalog (lAbdo et al.ll2010l) situ- 
ated within an angle < 10° from the blazar of interest. We model 
the broad band (0.1 GeV -10 TeV) spectra of the sources with 
a two-component model containing intrinsic emission from the 
primary source and emission from an electromagnetic cascade 
initiated by the absorption of very-high-energy gamma-ray in- 
teractions with the EBL, as explained in the next section. 



Name 


RA 


DEC 


z 


r 




lES 0229-H200 


38.203° 


20.288° 


0.14 


1.2 


5.0 


RGB J07 10+591 


107.625° 


59.139° 


0.13 


1.6 


1.0 


lES 1218+304 


185.341° 


30.177° 


0.18 


1.7 


2.5 



' Note that because of the difference in observation techniques, GeV 
and TeV measurements are accumulated at different time scales (year 
for the GeV band data and night-by-night in the TeV band). 



Table 1. Blazars considered in the analysis leading to the lower 
bound on EGMF. F and E^ut are the limiting values of photon 
index and the cut-off energy derived from the fit of the Fermi and 
TeV data with the direct emission (cut-off powerlaw attenuated 
by the pair production on the EBL) plus cascade model, under 
the requirement of a minimal cascade contribution. The cut-off 
energy in in units of TeV. 



For PKS 1424+240 the uncertainty of the source redshift 
does not allow firm predictions of the amount of cascade emis- 
sion to be made. Taking this into account, we exclude this source 



http :// http://fermi.gsfc.nasa.g0v/ssc/data/analysis/l 



A. M. Taylor, I. Vovk and A. Neronov: EGMF Constraints from Simultaneous GeV-TeV Observations of Blazars 



3 



from our analysis. For three sources, Mrk 501, 3C 66A, and 
PKS 2155-304, the slopes of the intrinsic spectra found in the 
fitting procedure are relatively soft, with photon indices F ^ 2. 
In such settings, the cascade emission gives a sub-dominant con- 
tribution to the source spectrum in the GeV band and no sensible 
constraints on the magnetic field along the line of sight could be 
derived from the analysis. The remaining three sources, listed in 
Table [1] have harder intrinsic spectra. We note that quite a small 
range of blazar redshifts is left, as seen in Table 1, resulting from 
the application of the selection criterias, mentioned above. For a 
given redshift bin size, there is more chance of finding a hard- 
spectrum blazar for the larger redshift values. On the other hand, 
for distant blazars absorption on the EBL in very significant and 
may cause them to become invisible in the TeV band. The com- 
bination of the above mentioned effects leads to the selection of 
the set of three sources, listed in Table 1 with quite a narrow 
redshift range. The analysis of the data on these three sources 
constrains the EGMF. 



3. Monte Carlo simulations 

Extragalactic electromagnetic cascades, in the presence of non- 
negligible EGMF (> 10"^" G), evolve both spatially and energet- 
ically as the propagation front of the emission moves away from 
the source. The nature of the constraints on the EGMF derived 
from the timing and imaging analysis of the signal produced by 
electromagnetic cascades developing in intergalactic space can 
be qualitatively understood from the decomposition of the cas- 
cade signal in space and time, which we illustrate using Monte- 
Carlo simulations. 

In order to convey the key features introduced by the spa- 
tial evolution of the cascade, we start with a consideration of 
simplified situation of a collimated primary gamma-ray beam 
(i.e. neglecting the finite jet opening angle). To further simplify 
the consideration, we consider an idealized situation in which 
a distant source injects primary gamma-rays at a fixed energy 
£() = 100 TeV. The energy £0 is chosen in such a way that 
the optical depth with respect to pair production on the EBL is 
t(Eo) » 1. The resulting spectra of the cascade emission de- 
pend only weakly on Eq as soon at t{E{)) » 1 . In these illustra- 
tive calculations we assume the source redshift z = 0.13, equal 
to the redshift of RGB J0710H-591. Later in this section we in- 
troduce non-zero intrinsic jet opening angle Oj,,, > and broad 
band emission spectrum and describe how our simplified results 
are altered. 

Fig. 1 shows the results of the model calculation. The arriv- 
ing spectra of the cascade emission, for the case in which the 
cascade develops in the presence of a negligible EGMF, is the 
long dash, short dash line seen in both panels of figure [1] This 
arriving photon spectrum may, in fact, be obtained simply us- 
ing a "kinetic equation" description. Indeed, we have compared 
our results using a Monte Carlo descrip t ion us ed in this work 
(for more details see lTavlor & Aharonianl(l2008h ) and the kinetic 
description results for t he case of a negligible B-field, used in 
iNeronov & Vovkl (l201d) . finding good agreement in all cases. 

For these results and the results th roughout this paper, the 
EBL model of lFranceschini et al] (l2008h is adopted, in which the 
evolution of the EBL with redshift is accounted for Furthermore, 
the cosmological parameters of Hq - 70 km s"' Mpc"', - 
0.3, and Qa - 0.7 are also adopted. 

In the calculations throughout this work, the EGMF are as- 
sumed to be describable by patches of coherent (uniform) mag- 



RGB J0710+591 (z=0.13) 
6=3x10"'"' G 




Fig. 1. The arriving energy fluxes following the injection of a 
lO'** eV photon flux from a source at redshift z - 0.13, with 
an intervening EGMF=3 x 10 '^G. Top: The arriving flux is de- 
composed into that observed from different observer angles, 0obs- 
Bottom: The arriving flux is decomposed into that arriving with 
different delay times, fdeiay In both plots, the long-short dash 
line represents the envelope flux containing the non-decomposed 
spectrum. 
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Fig. 2. The mean time delay incurred by a photon energy flux fol- 
lowing its injection at 10'"* eV from a source at redshift z = 0.13, 
for different intervening EGMF, in the range 10"^-10 G. 



netic field, each patch being Ab - 1 Mpc in sizqj. Different 
magnetic patches have their fields orientations chosen com- 
pletely independently (randomly) of each other. Such a descrip- 



^ We have verified that the derived limits on EGMF do not depend 
on /Ig for /Ig > 1 Mpc by comparing the results of calculation of the 
cascade for Ag = I Mpc with those for Ag = 30 Mpc. 
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Fig. 3. A diagram depicting the flaring of the electromagnetic 
cascade development in the presence of extragalactic magnetic 
fields. The initial (conical) jet emission whose power feeds the 
electromagnetic cascade is also shown, with the inner shaded 
region on the right representing the region filled by the intrinsic 
cone emission. 



tion amounts to assuming that the magnetic field power spectrum 
places a significant proportion of the magnetic field energy den- 
sity at the longest length scales. Since, if this is not the case, the 
magnetic field deflections will be somewhat weaker, our results 
obtained in this work under this assumption can be considered 
conservative. 

The presence of a non-negligible magnetic field (> 10"^° G) 
introduces both spatial spreading of the cascade front in angle, 
away from the initial beam direction, and the sig nificant grow th 
of the cascade front depth (arrival time spread (lPlagal[l995l) ). 
Spatial and time decomposition of the full cascade spectrum 
arising at non-zero EGMF is shown in the top and bottom panels 
of Fig.[Tl respectively. 

For the calculation shown in the top panel of Fig. [T] an 
observer is placed at an off-axis angle 9obs with respect to the 
gamma-ray beam, at a distance D from the source. The observer 
detects only photons incident on the sphere of the radius R - D 
at small incidence angles 6 < Opsp with respect to the normal to 
the sphere, flpsp corresponds to the point spread function of ob- 
server's telescope. An approximation of Opsp for the FermijhP^ 
instrument, valid for photon energies in the range 30 MeV to 
300 GeV 



0PSF ~ 1-7° (fiy.Gev) 



-0.74 



1 



15 



0.37 



(1) 



where fiy.cev is the photon energy in GeV, is used for the pur- 
poses of our calculation^ Bin widths of 0.1° (i.e. photons are 
collected from a ring 9obs ±0.1°) are used for the observer posi- 
tion. In the bottom panel of Fig. [1] only photons arriving with a 
particular delay time, relative to that of the propagation time of 
the direct gamma-ray photons from the source, are shown. The 
width of the time bins is half a decade around the reference fdeiay 
value. 

Fig.[2]shows the time delay of the cascade photons as a func- 
tion of the photon energy for different magnetic field strengths. 
Coarse bin widths of 1 per logarithmic decade in time are used 
for time delay binning. One can see that the time delay de- 
creases proportionally ~ or ~ E^^, with the increasing 
photon energy Ey. in agreement with the analytical calculations 
(Neronov & Semikoz 2009). Change of the the slope of the en- 
ergy dependence of the time delay at high energies is related to 



the energy dependence of the mean free path of the primary y- 
rays. 

3.1. Jet Opening Angle 

A realistic description of jet cascade effects in an EGMF should 
take into account the intrinsic spread in angles of the injected 
gamma-rays, within the opening angle of the blazar jet, fljet as 
well as the spread in initial energies of gamma-rays, £0. 

To take into account the finite opening a ngle of the b l azar je t 
we use the same procedure as described in lElviv et all (|2009|) . 
utilizing zero jet width results. Using a random remapping of 
the arriving particles' position and velocity vectors on the arrival 
surface, of the zero jet width results, the effect of the jet width on 
the cascade results is obtained. Fig.|3]illustrates the geometrical 
set up of our Monte Carlo calculations. In what follows we fix 
the jet width to 0jet - 5°, which is typical opening angle of the 
jets observed in radio galaxies and is the angle fljet - F^r' inferred 
from the typical bulk Lorentz factors of the blazar jets Fjet ~ 10 
derived from the gamma-ray observations. We (arbitrarily) place 
the observer at Gobs = 2°. 

Cascade spectra for the case of a non-zero jet opening angle 
can be qualitatively understood as the weighted sums of different 
components of the decomposition of the full cascade spectrum 
shown in Fig. [1] (i.e. adding cascade emission produced by pri- 
mary photons emitted in different directions within fljet is equiv- 
alent to displacing the observer to a different 0ob,s with respect 
to the primary beam direction). The effect of introducing an in- 
trinsic jet opening angle is, therefore, to wash away the small 
angle/small time delay information present in the high energy 
emission. 



3.2. Intrinsic spectra of the blazars 

It is conventionally assumed that the spectra of the blazars mea- 
sured in the Fermi energy band are "intrinsic" spectra, produced 
at the source and that only the spectrum above Ey>\QQ GeV is 
strongly influenced by the effects of pair production on the EBL. 
This assumption does not work, a priori, if the uncertainty of 
the strength of EGMF is taken into account. In particular, if the 
EGMF is close to zero, cascade emission can give a significant 
contribution to the observed source spectrum in the GeV energy 
band, if the intrinsic spectrum of the source is hard, so that most 
of the power is initially emitted at energies Ey > 100 GeV. In 
particular, it is even possible that the cascade component gives a 
dominant contribution in the source flux in the GeV-TeV band, 
as was demonstrated for the pa rticular case of the blazar Mrk 
501 by Aharonianet al] (12001 1) . More generally, the observed 
spectrum in the 0.1 GeV - 10 TeV energy range, simultane- 
ously measured by Fermi and by the ground-based Cherenkov 
telescopes is a superposition of the direct source flux (with an 
a-priori unknown spectrum) and a cascade flux. The spectrum 
of the cascade flux can be unambiguously derived from the in- 
trinsic source spectrum, if a particular configuration of EGMF 
is assumed. Thus, proper modeling of the broad band y-ray data 
of blazars in the 0.1 GeV - 10 TeV energy range has to be used 
to derive both the intrinsic source spectra and the properties of 
EGMF from the model fits to the observed spectra. 

We limit the choice of models for the intrinsic source spectra 
with the cut-off powerlaw type models 



See h ttp://www-glast.slac.stanford.edu/software/IS/| 
glast_lat_performance.htm 



dNy _^ I Ey\ 



(2) 
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Fig. 4. The arriving photon energy flux, resulting from the "minimal cascade" model for Mrk 501 (top), PKS 2155-304 (middle) 
and 3C 66A (bottom). Long-dash-dotted line shows the intrinsic source spectra. Red solid lines show the spectra with account of 
the cascade component. 
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Fig. 5. Top: The arriving photon energy flux, resulting from the "minimal cascade" model for RGB J0710+591 (left), lES 0229+200 
(center) and lES 1218+304 (right). Fermi upper bounds are at the 99% confidence level. Notations are the same as in Fig.|4| Blue 
dash-dotted line shows the cascade component of the spectrum. Bottom: the results for the "maximal cascade" model. 



described by two para meters, the photon index F and the cut-off 
energy Scut- Following lNeronov & Vovl3 ( l2010l) . the values of F 
and ^cut are chosen such that the absorbed spectrum, along with 
the subsequent cascade contribution, give both good fits to the 
complete multi-wavelength (GeV and TeV) data set and mini- 
mize the cascade contribution. Technically, if the cut-off energy 
iicut is not constrained by the data, the lower bound on ^cut is de- 
rived from the TeV data and this lower bound is assumed for the 
derivation of the cascade component of the spectrum. Since this 
procedure minimizes the fraction of TeV flux absorbed during 
propagation, the fit obtained is labelled as our "minimal case" 
model for the cascade contribution. 

In addition to the powerlaw type spectra we also consider 
a qualitatively different model for the intrinsic source spectra. 



which results in the dominant cascade contribution in the GeV 
band. We call this case the "maximal" cascade model. As is men- 
tioned above, as soon as the bulk of the initial source power 
is injected at energies far above the threshold for pair produc- 
tion on the EBL, the resulting shape of the cascade compo- 
nent of the spectrum almost does not depend on the details of 
the intrinsic source spectrum. Taking this into account, we ar- 
bitrarily fix the intrinsic photon index to F = 1 and the cut-off 
energy to Ecut — 100 TeV. The results for the observed GeV- 
TeV band spectr um in the maxim a l casc ade model are similar to 
those derived bv I Aharonian et al" I (l200l for the case of negligi- 
ble EGMF. 
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4. Results 

4.1. Minimal and maximal cascade model fits to the spectra 

From the combined broad 0.1 GeV - 10 TeV band gamma-ray 
spectra of the sources listed in Table [T] we find that in general, 
the high energy end of the Fenni/LAT spectra at Ey ~ IQQ GeV 
perfectly matches the low energy end of the spectra obtained 
with the ground based gamma-ray telescopes. This indicates that 
no additional "intercalibration factors" between different instru- 
ments are needed in the modeling of the source spectra. 

The observed 0.1-100 GeV band spectra of Mrk 501, 
PKS 2155-304 and 3C 66A have photon indices close to T = 2 
with almost no curvature of the spectrum over the entire energy 
range. Such a "flat" spectral energy distribution is difficult to ob- 
tain in the "maximal cascade" model. Although the "universal" 
cascade spectrum has average photon index close to F = 2 in 
the 0.1-100 GeV band, its deviation from the pure powerlaw be- 
havior is significant, so that this model is in contradiction with 
the observed spectrum for these blazars. This suggests that the 
dominant contribution to the observed source spectrum in the 
cases of Mrk 501, PKS 2155-304 and 3C 66A comes from di- 
rect gamma-rays from the source as shown in Fig. |4] 

Addition of the cascade contribution to the intrinsic emis- 
sion component in Mrk 501, PKS 2155-304 and 3C 66A leads 
only to a moderate modification of the source spectrum in the 
GeV energy range. Slight adjustment of the intrinsic spectrum 
photon index is sufficient to make the observed spectrum con- 
sistent with the model even for the case of zero EGMF strength. 
The reason for the sub-dominance of the cascade contribution 
in these sources is the softness of the intrinsic source spectrum 
(F ^ 2). The fraction of the source power absorbed on the way 
from the source to Earth is at most comparable to the primary 
source power at energies Ey < 100 GeV, so that cascade emis- 
sion can not dominate over the intrinsic source power. 

The situation is different, however, for the case of 
RGB J0710-H591, lES 0229H-200, and lES 1218-H304. The in- 
trinsic spectra of these sources are harder than F = 2. The hard 
intrinsic spectrum of the source leads to a large power output at 
energies above 100 GeV, even in the case when a possible high 
energy cut-off in the source spectrum is taken into account. In 
the "minimal" model shown in Fig.|5] the intrinsic source power 
output in the TeV energy band is an order of magnitude higher 
than that in the GeV band. A significant fraction of the of the 
intrinsic source power in the TeV band is absorbed through pair 
production on the EBL. This absorbed power is re-emitted in the 
GeV band, so that the cascade emission in the 0.1-10 GeV range 
dominates over the intrinsic source emission (solid red curve in 
FigE). 

The fit of the model (|2]i alone to the observed spectra re- 
sult in reduced close to 1 both for RGB J07 10+591 and 
lES 1218+304. For lES 0229+200 only an upper bound on 
the source flux in the GeV band could be derived from Fermi 
data. Account of the cascade contribution to the source spec- 
tra, which has to be present for the case of zero EGMF, vi- 
olates the upper bound on the GeV band source flux in the 
case of lES 0229+200. In the case of RGB J0710+591 and 
lES 1218+304 account of the cascade contribution leads to the 
worsening of the fit to the spectrum from = 8 (7 d.o.f.) to 
= 18 (RGB J0710+591) and from ^ ^ (jg (J ^ q 
= 57 (lES 1218+304). This implies that the model with 
cascade component calculated under the assumption of zero 
EGMF is ruled out at 98.8% and > 99.99% confidence levels 
for RGB J07 10+591 and lES 1218+304, respectively. 



RGBJ0710+591 (z=0.13) Injection spect|ym 
minimal case ^^-ir^ 




RGB J0710+591 (z=0.13) Injection spectrum 

I in 




RGBJ0710+591 (z=0.13) Injection spectrum 




8 9 10 11 12 13 14 



loQioEy [eV] 

Fig. 6. The effects of the presence of an EGMF on the arriving 
cascade from the blazar RGB J07 10+591. Two top panels show 
the results for the case of the "minimal cascade" model of the 
y-ray spectrum, for the two possibilities of suppression of the 
cascade flux via extension of the cascade source and time delay 
of the cascade emission. Two bottom panics show the results for 
the "maximal cascade" model. 



In the case of lES 0229+200, the minimal model parame- 
ters suffer from a large uncertainty: both the cut-off energy ^cut 
and the s pectral slope F could no t be derived from the data. 
Following iNeronov & Vovkl (120101) we fix the model parame- 
ters for this source in such a way that the total flux of the cas- 
cade component of the source spectrum is minimized. This is 
achieved with a very hard value of F ~ 1.2 and relatively high 
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1ES 0229+200 (z=0.14) Injection spectrum 
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Fig. 7. Same as in Fig.|6]but for lES 0229-3+200. 



Fig. 8. Same as in Fig.|6]but for lES 1218+304. 



cut-off energy Ecut ~ 5 TeV, so that the "minimal" and "maxi- 
mal" models for lES 0229+200 are not very different. 

The problem of inconsistency of the predicted cascade flux 
with Fermi measurements in the GeV band is encountered in the 
"maximal" model shown in Fig. |5] Assuming that the cascade 
emission dominates over the direct source emission in the TeV 
band, one finds that the expected cascade flux level in the 0.1- 
10 GeV band is higher than the observed source flux, if EGMF 
strength is B = (solid red curves in Fig.|5]l. 

4.2. Implications for EGMF 

The fits to the combined GeV- TeV band spectra found in the 
minimal and maximal models for the cascade contribution under 



B = assumption are shown in Fig. |5] for RGB J0710+591, 
lES 0229+200, and lES 1218+304. For these sources, both 
minimal and maximal models, calculated under the assumption 
of zero EGMF, are not acceptable. The main source of discrep- 
ancy between the model predictions and the data is an over pre- 
diction of low energy flux, Ey < I GeV, due to the presence of 
the cascade contribution to the total source flux. 

As discussed in section[3] a decrease in the cascade contribu- 
tion to the arriving flux at low energies can be achieved through 
the introduction of non-zero magnetic field in the cascade devel- 
opment region. The effect of a non-zero magnetic field on the 
cascade component is shown in Figs. |6] 111 and[8]for the cases of 
suppression of the cascade flux due to the large extension of the 
cascade source and due to the time delay of the cascade signal. 
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Fig. 9. Bounds on magnetic field derived from the simultaneous 
GeV-TeV data. Blue shaded regions show the previously known 
bound s on B and Ab, summarized by Neronov & Semikoz: 
(I2009h . Orange shading shows the upper bound on B,Ab which 
co uld be generated before the e poch of recombination, derived 
bv lBaneriee & Jedamzild ( l2004l) . 



For the case of suppression due to extended nature of the 
cascade source, the presence of magnetic field modifies the 
cascade spectrum at GeV energies only if the magnetic field 
strength is B > 10 G, for the case of RGB J07 10+591, 
lES 0229+200, and lES 1218+304. The minimal magnetic field 
strengths needed to make the model source spectra consistent 
with the data can vary between lO"'*" and 10"'^ G, depend- 
ing on the adopted model source (from "minimal" to "maxi- 
mal", through all the "intermediate" possibilities) and the model 
of the EBL. The tighest bound is derived from the data on 
IBS 0229+200, at the level of 10"'^ G, which is consistent 
with the bounds found under simil ar assumptions about th e 
cascade suppressi o n mechanism by iNeronov & Vovkl (l20Tol) : 
iDolag et al.1 OOld) : iTavecchid (l2010l) . We stress that the bound 
should be considered only as an order-of-magnitude estimate, 
due to the significant uncertainty of the shape and overall nor- 
malization of the cascade introduced by the uncertainty of the 
normalization and spectral shape of the EBL. 

For the case of suppression of cascade emission due to the 
time delay of the cascade signal, one assumes that the primary 
source is active only during a limited period of time, just about 
the time span of gamma-ray observations (fsource ~ 1 yr). Time 
delay of the cascade signal by fdeiay > ^source would lead to the 
suppression of the cascade flux by a factor fsouice/ fdeiay Figs.|6l 
|7J and [8] show that time delay of the cascade signal starts to in- 
fluence the cascade emission signal at GeV energies when the 
magnetic field strength reaches ~ 10""^ G. Similar to the case of 
suppression due to the extended emission, the precise value of 
B necessary to suppress the cascade emission depends on the 
adopted source and EBL models. The tighest lower bound is 
again derived from the data on lES 0229+200, at the level of 
10 G. This bound should also be considered as an order-of- 



magnitude estimate because of the remaining uncertainty in the 
measurements of the spectrum of EBL. 

Note that the bound B > 10"'^ G derived assuming sup- 
pression of cascade emission due to the time delay of the cas- 
cade signal in the case of lES 0229+200 is by 1.5 orders of 
magnitude stronger than t he bound derived fro m a similar anal- 
ysis of the same source bv lDermer et al We believe that 
the main source of discrepan cy between the result obtained in 
the present work and that of iDermer et al.l (1201 Ol) is the sim- 
plified analyt i cal tre atment of the cascade emission adopted by 
iDermer et al.l ( l2010h . The simplified treatment of the cascade 
has led to an under-estimate of the cascade flux at high ener- 
gies Ey > 10 GeV and an over-estimate of the strength of sup- 
pression of the cascade emission due to the time delay at low 
energies Ey < 10 GeV. 

Furthermore, we note that our Hmit of B > 10 G from 
the time delay of the cas cade signal is consistent with the results 
of a similar analysis by iDolag et"an (1201 Oh . who found some- 
what tighter bound B > 10 G, assuming a larger minimal 
possible time delay, fjeiay > 100 yr in the cascade emission from 
lES 0229+200. 

A summary of the Umits on magnetic fields in the inter- 
galactic medium, which can be derived from the simultaneous 
GeV-TeV band observations is shown in Fig. |9l In our analysis 
we have considered the bound on the EGMF strength assum- 
ing a fixed magnetic field correlation length Ab - I Mpc. If 
the EGMF correlation length is Ab ^ 1 Mpc, the lower bound 
on EGMF strength does not depend on Ab because the cooling 
distance of e^e" pairs is much shorter than the typical size of 
regions in which EGMF is correlated. We have explicitly ver- 
ified this by making a control run of Monte-Carlo simulations 
with Ab = 30 Mpc and comparing the results with the case 
Ab - I Mpc shown above. On the other hand, if Ab < I Mpc, the 
inverse Compton cooling distance becomes larger than the size 
of the regions with correlated EGMF. This means that electrons 
and positrons pass through regions with different magnetic field 
orientations during their cooling. As a result, the deflection an- 
gle scales proportionally to the square root, rather than linearly 
with the propagation distance on the distance scales comparable 
to the inverse Compton cooling length. This explains the im- 

— 1/2 

provement of the lower bound on the EGMF strength B ~ A^ ' 
at Ab < 1 Mpc; stronger magnetic field is required to deviate 
electron trajectories by a given angle. 

5. Conclusion 

In this paper we have derived constraints on the strength of mag- 
netic fields in the intergalactic medium from simultaneous obser- 
vations of blazars in the GeV band (by Ferm//LAT telescope) and 
TeV band (by ground-based y-ray telescopes). The constraints 
stem from the requirement that the GeV band signal from elec- 
tromagnetic cascade initiated by the absorption of the primary 
TeV y-rays in interactions with Extragalactic Background Light 
should be suppressed by deflections of electron-positron pairs by 
magnetic fields in the intergalactic medium. Non-observation of 
the cascade emission by Fenni/LAI telescopes imposes a lower 
bound on the cascade flux suppression factor which could be 
converted to a correlation length dependent lower bound on the 
strength of magnetic field. 

We have found that constraints on the magnetic field 
strength could be derived from the y-ray data on three blazars, 
lES 0229+200, RGB J0710+591 and lES 1218+304 (out of 
seven, for which simultaneous GeV- TeV data are available). For 
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all three sources, we have performed detailed modeling of the 
spectral characteristics in the broad (0.1 GeV to 10 TeV) en- 
ergy range. We have fitted the observed y-ray spectra with a 
two-component spectral model which consists of both direct ab- 
sorbed emission from the blazar and a cascade emission compo- 
nent calculated using detailed Monte-Carlo simulations of the 
cascade development. The observational data are inconsistent 
with the models in which the cascade emission is calculated as- 
suming zero magnetic field strength in the cascade development 
region (extending to ~ 100 Mpc distance from the source along 
the line of sight). 

The minimal magnetic field strength required to achieve suf- 
ficient suppression of the cascade signal depends on the assump- 
tion about the mechanism of suppression of the cascade signal. 
If the suppression is due to the time delay of the cascade emis- 
sion following a period of enhanced activity of the source in the 
TeV band (with duration ~ 1 yr), then the minimal required field 
strength is B ~ 10 G in the case of the field with large cor- 
relation length Ab > I Mpc. If the (unknown) correlation length 
is Ab < 1 Mpc, the minimal needed magnetic field strength is 
larger by a factor (Ab/1 Mpc)"'^^. If the suppression of the cas- 
cade emission is due to the fact that the size of the cascade source 
is much larger than the point spread function of Fermi/LAI tele- 
scope, rather than due to the time delay of the cascade signal, the 
lower bound is B > 10 G, with the same dependence of the 
correlation length Ab- 

The two possibilities for suppression of the cascade emission 
could be distinguished via a search of the delayed GeV y-ray 
emission following strong TeV band flares of blazars or of the 
extended emission around TeV blazars in the GeV- TeV band. 
The distinguishing feature of the delayed cascade emission is 
the characteristic energy dependence fdeiay ~ £y~'^ or ~ is^-, as 
shown in Fig. |2] If the real magnetic field strength is close to 
B ~ 10 G, the extended emission around TeV blazars should 
be detectable with Fermi, while the time delay of the cascade 
emission might be detectable at higher energies (~ 100 GeV) by 
the ground-based y-ray telescopes. 

We have investigated the dependence of the derived limit on 
the assumptions about the intrinsic spectrum of the sources by 
considering the extreme cases of "minimal" and "maximal" cas- 
cade contributions. In the "minimal" case the parameters of the 
intrinsic spectrum of the source are chosen in such a way that 
they minimize the total flux in the cascade component (at zero 
magnetic field strength). In the "maximal" cascade model the 
cascade flux dominates over the intrinsic source flux at all ener- 
gies up to the TeV range. Surprisingly, the lower bound on the 
magnetic field strength is practically independent on the choice 
of the model used to fit the observed GeV - TeV band spectra. 
This is explained by the fact that the amount of power trans- 
ferred by the cascade from the TeV to the GeV energy band is 
determined only by the measured TeV flux from the direction of 
the source and is not sensitive to the origin of the TeV y-rays (if 
they are intrinsic to the source or produced in the course of de- 
velopment of electromagnetic cascade close to the source). Since 
"minimal" and "maximal" cascade models represent the two ex- 
treme possibilities for the possible amount of cascade contribu- 
tion to the source flux, we conclude that the uncertainty of the 
modeling of the observed source spectrum introduces an uncer- 
tainty by a factor of ~ 1 in the derived lower bound on EGMF 
strength. 
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